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Abstract 



^s^ ' We examine the parameter accuracy that can be achieved by advanced 

^ , ground-based detectors for binary inspiralhng black holes and neutron 

fT^ ■ stars. We use the 2.5 PN spinning waveforms of Arun et al. [10]. Our 

("^ ' main result is that the errors are noticeably different from earlier studies. 

^O , An important contribution to this difference comes from self-spin terms 

>0 ' at 2PN order not previously considered. While the masses can be deter- 

fT^ I mined more accurately, the individual spins are measured less accurately 

("^ . compared to previous work, in some cases by more than a factor of ten. 

O^ ' We also examine several regions of parameter space relevant to expected 

sources and the impact of simple priors. A combination of the spins is 
measurable to higher accuracy and we examine what this can tell us about 
spinning systems. 

X' 

.•?.■ 1 Introduction 

Upgrades to gravitational wave interferometric detectors are currently being 
performed [1]. Once completed these advanced detectors are expected to de- 
tect tens of astrophysical signals per year [2] . The inspiral of massive compact 
bodies in binary systems is one of the most important expected signals. Detec- 
tion of the gravitational waves from such events can provide information about 
the parameters of the compact bodies involved. Knowing the masses of the 
objects involved is central to classifying the observed objects as neutron star 
or black hole candidates. In addition, the objects may have intrinsic spin and 
General Relativity predicts that this will affect the motion of the bodies and 
the gravitational waveforms produced at a measurable level [3] . 



To detect the signals in the background noise the method of matched fil- 
tering is employed [4]. Theoretical templates based on Post-Newtonian (PN) 
expansions in General Relativity are employed to search the detector output for 
signals. The goodness of fit of a particular template to the data can be most 
easily seen by computing a signal-to- noise ratio (SNR). Because the detector 
is noisy, the best fitting set of parameters may not be the same as the true 
parameters of the source. The error in the estimation of the source parameters 
depends in part on the signal to noise ratio. A theoretical indication of the pa- 
rameter errors can be computed by simulating signals and multiple realisations 
of the noise or by using the covariance matrix of Fisher information theory. We 
employ the Fisher information method here. 

As the true waveform in General Relativity is approximated by a Post- 
Newtonian expansion, a systematic error is introduced when truncating the 
true waveform. From a theoretical viewpoint, it would be advantageous if this 
systematic error is smaller than other sources of noise. The order to which the 
Post-Newtonian expansion is needed in order to provide a satisfactory model for 
detection and parameter extraction is still unknown, although there are some 
hints of at least a lower bound [5] [6] . 

At the moment we do not possess PN expansions beyond the level that 
seem satisfactory and so it is prudent to use the highest order PN expansion 
available. Previous work has examined the measurement sensitivities of ground 
based interferometers for 1.5PN [7] and 2PN [8] [9] orders. Here we extend those 
works to the recently available 2.5PN spinning waveforms of Arun et al. [10] 
which also includes self-spin terms at 2PN not considered in previous work. Non- 
spinning waveforms are available to 3.5PN order and their parameter estimation 
has been examined in [11]. 

As noted in [13] the gravitational wave flux for test particle non-spinning 
inspirals is badly behaved at 2.5PN level. Table I of [11] gives parameter esti- 
mation for non-spinning binaries at increasing steps of 0.5 in the PN expansion 
up to 3.5PN order. Although the convergence of the series looks very slow, the 
parameter errors for non-spinning waveforms shows no large fiuctuations at the 
2.5PN order. 

We concentrate here on aligned spinning systems. With the assumption of 
aligned spin and orbital angular momentum axes, the waveforms depend on six 
parameters. Two of the parameters, tc and <f>c^ are extrinsic parameters deter- 
mining the time and phase of the wave at coalescence. The point of coalescence, 
which is not well modeled by the PN expansion, could be replaced by any other 
chosen point in the inspiral - both parameters are effectively integration con- 
stants. Although the template waveform match does depend on their value, the 
parameter error ranges do not. 

The other four intrinsic parameters are the masses and spins of the two 
inspiralling objects. It is ultimately only the masses and spins of the objects 
that we are interested in determining, but these parameters are correlated with 
tc and (j)c- We will examine the effect of redefining the coalescence phase, (/)c, 
and also the effect of imposing 27r periodic priors on this parameter. 

Three important differences exist between the 2.5PN waveforms investigated 



here and earlier work for ground based detectors. Firstly the extension to 2.5PN 
order provides seven expansion coefheients in the waveform expanded as a series 
in the frequency, with only six parameters to be measured. These expansion 
coefficients are sometimes called "chirp times" in the time domain waveforms. 
Previous work at 1.5 FN order for ground based detectors [7] considered five FN 
expansion coefficients depending on five parameters: ic, (f'c, the reduced mass, 
/It, chirp mass, A4 and /3, a parameter related to the spin. Frevious work at 
2FN order [8] considered six expansion coefficients with six parameters - adding 
an extra spin related parameter, a. Adding more parameters typically reduces 
the accuracy at which any one parameter can be measured, but adding more 
expansion coefficients can counteract this effect. At 2.5FN order we are adding 
another coefficient without increasing the number of parameters, since the extra 
coefficient can be expressed as a combination of the others. 

Secondly, the introduction of new terms independent of the frequency in the 
waveform phase greatly increases the correlation between the phase at coales- 
cence, (j)c and the other parameters, in particular the spins. Errors in (j)c can 
then greatly magnify errors in the spins. This makes it even more imperative 
that errors in the tpc are controlled in order not to over-contaminate the errors 
in the mass and spin parameters. 

Thirdly, the waveform considered here contains self-spin interaction terms of 
the type spin(l)spin(l) and spin(2)spin(2) at 2FN order, that were not consid- 
ered in previous analyses. The 2.5FN waveform of [10] contains these self-spin 
interaction terms at 2FN order, which are not present in the 2FN waveform used 
in [8]. These terms were first calculated in [12]. A major part of the increase in 
spin uncertainty can be traced to these self-spin interaction terms. The 2.5FN 
waveform used here, only contains spin-orbit terms at 2.5FN and it remains to 
be seen whether further, as yet uncalculated FN terms, will affect our results. 

The first issue to be addressed in our work is; do the parameter errors dif- 
fer significantly from results obtained using lower order FN expansions? The 
answer to this question will provide some guide to the question of what or- 
der of FN expansion is sufficient to model the true General Relativity signal 
to an acceptable level. A second question we will examine is; at this level of 
FN approximation and with the assumption of aligned spin and orbital angular 
moment, does intrinsic spin matter? To what extent can the experiments dis- 
tinguish a spinning system from one without spin? A further, slightly different 
question is; how accurately can the expected data constrain the individual spins 
of each of the inspiralling objects? 

For spinning black holes whose exterior is exactly described by the Kerr 
metric a theoretical upper limit M^ > J exists bounding the magnitude of spin 
angular momentum, J, to be less than the square of the mass, M , in natural 
geometric units. Beyond this limit the object can no longer be described as a 
Kerr black hole but a naked singularity. This relationship can be expressed in 
terms of the spin parameter x = J/M^ as x ^ 1- To truly claim that a massive 
object is a black hole one must show that the spin of the object satisfies this 
bound. 

Black holes can be born with intrinsic spin or be spun up due to the accretion 



of matter. Theoretical studies indicate that the maximum amount of spin that 
a black hole can accumulate due to accretion corresponds to x = 0.998 [14]. A 
number of highly spinning black holes, with spin parameters x > 0.9 may have 
been found using X-ray techniques [15]. 

Neutron stars are expected to have maximum spins beyond which they break 
up due to centrifugal forces [16]. The x = 1 case for a neutron star corresponds 
to a frequency in the lO'^Hz range. Ultra-fast millisecond pulsars have been 
observed with spins up to 700Hz. The three known binary pulsar systems that 
will coalesce in a Hubble time all have spins corresponding to x < 0.02 [17]. 
Non-compact objects, such as the Earth or Sun, actually do violate the x = 1 
bound, but it is not expected that such non-compact objects will be able to 
survive the tidal forces present in the last moments of inspiral in the frequency 
ranges of interest to us. 

We will attempt a wider investigation of the 2.5 PN order parameter space 
by allowing the inspiralling objects to be spinning and examine the effect that 
this has on the parameter estimation. The errors in the parameters depend 
sensitively on the region of parameter space that is probed. Covering the full 
parameter space in detail is time-consuming and leads to the problem of pre- 
senting the resulting twelve dimensional space - the six parameter errors for the 
six parameters. We will instead pick a few areas of the full parameter space that 
seem of particular interest and present the results mainly in a set of tables for 
ease of comparison. The Fisher matrix formalism, if reliably implemented, is far 
more efficient at probing large parts of the parameter space than full Markov 
chain Monte Carlo simulations. 

We also examine the effect of allowing the cndpoint of inspiral to change 
for spinning objects as has been suggested by Hanna et al. [18]. Previous 
studies have conservatively adopted the Schwarzschild IS CO as the endpoint of 
the inspiral. The ISCO frequency of an extremal Kerr black hole is roughly a 
factor seven larger than the ISCO frequency of a Schwarzschild black hole with 
the same mass. This may allow us to extend the range of validity of the PN 
waveform and we investigate what effect this has on the parameter estimation. 

We focus on determining the intrinsic parameters of mass and spin. Param- 
eters of the source objects such as the distance to the source and the orientation 
of the binary will affect primarily the amplitude and we will ignore these affects 
here, effectively by normalising the amplitude to obtain a given signal-to-noise 
ratio (SNR) and assuming the binary is directly above the interferometer, with- 
out any Doppler shift relative to the detector. Here we normalise the SNR to be 
10 for a single detector. This is slightly above the threshold for detection and 
compatible with previous methodology. This choice is merely a normalization 
of the gravitational wave amplitude and results for other SNR values can be 
found by a simple linear scaling. The SNR for which the Fisher matrix formal- 
ism agrees with Monte Carlo simulations may lie above 10 [19], but this is the 
conventional scale at which to compare results. Estimates for the performance 
of multiple detectors operating together can be indicated by combining the in- 
dividual SNRs in quadrature [7]. Multiple detectors have added advantages in 
dealing with realistic non-Gaussian non-stationary noise, which we do not treat 



here. 

If the compact objects have spin axes that are not ahgned with the orbital 
angular momentum the axes will precess leading to a modulation of the am- 
plitude and other effects [3]. For simplicity we will also ignore this here. 2PN 
waveforms, without spinlspinl self-spin coupling terms, have been investigated 
for the LISA space detector for the fully aligned case in [20], precessing with 
random spin orientations in [21] and for the partially aligned case in [22]. Fur- 
ther parameters of the source such as mass and energy accretion rates, mass 
multipole distributions and internal equations of state will also affect the gravi- 
tational waveform. The mass multipole correction for spinning neutron stars is 
second order in the Post- Newtonian approximation [7], but we leave a further 
investigation of these effects to future work. 

2 Measurement error 

We follow the standard Fisher information methodology of [7] and [8] . Wc focus 
on the case of a single detector, with stationary, non-Doppler shifted, overhead 
source, for simplicity. With the assumption that the noise in the detector is 
stationary Gaussian with correlation function Cnir), for time separation r, the 
probability that the data takes some particular, discretely sampled, pure noise 
value. Si, given an assumption that no signal is present is 

In the continuum limit, using the Wiener- Kinchin theorem to relate the one- 
sided power spectral density, S'„(/), to the correlation functions, Cn{T), this 
can be written as 

In practice we will assume that the noise is dominant below some lower frequency 
bound /. and we will not follow the inspiral waveform beyond an upper frequency 
/, , so we can define an effective inner product such that 

J f. ^n[J) 

Now we can ask what is the probability of inferring a set of parameters, 9, giving 
rise to a signal, h{9), from a measured detector output, s, if there is a true signal 
ho in the data, such that s ~ ha + n, with n the noise. To do this we use the 
relation p{s\h{9)) = p{s — h{6)\0) and ask what is the probability of obtaining 
an output s — h{9), if only noise is present. A generic gravitational waveform 
strain amplitude, h{9), can be expanded around a given reference waveform, Hq, 

by 



where A6i denotes the difference between the ith parameter and the value it 
takes in Hq. Then the hkehhood, p{s\9), is given in stationary Gaussian noise 

by 

pis\e) ^ ^-{s-h(e),s-h(e))/2 

^ g(-(n,ri)/2+Aej(nJij)-AejAefc(ftj,ftfc)/2+...) /gx 

with the assumption that the data, s, is given hy s — ho + n and hi is defined as 
hi — dh/dO'^. In the Linear Signal Approximation (LSA) all higher derivatives 
of h can be dropped, and thus terms in the ellipses can be neglected. The term 
with (n, n) is just an overall normalisation and for noise with zero mean the 
term (n, hj) vanishes, so applying the LSA, corresponding to the high Signal to 
Noise Ratio (SNR) [23], the likelihood is just 

p{s\e) - e-''^'''"('^^-^'"=)/2, (6) 

and the posterior probability density function for the parameters 9^ can be found 
from p{9)p{s\9). The expectation of the (co)variance of the parameters, in an 
n parameter model, can be found from the inverse of the Fisher-information 

_ j M,Mjp{9)p{s\9)d" 
J p{9)p{s\9)d'^9 

where the angled brackets on the left hand side denote expectation. In the limit 
of trivial priors p{9), that are flat over the entire real line, this becomes 



(A^.A^,) . '^T:r;::z::r\ (7) 



{A9,A9j)^{K,h,)-\ (8) 

The Fisher information matrix is generated from the given PN waveforms using 
the above inner product by 



P . dh 

J- 1 



- f-J 



Si- (^) 



There are several ways of including the effect of priors on the posterior distribu- 
tion. Most easily treated are the Gaussian priors used by Cutler and Flanagan 
[7] and Poisson and Will [8]. These can be incorporated simply by adding a 
prior matrix to the Fisher matrix with terms given by l/ug, with ag the desired 
prior variance on the parameter 9. In this case the Gaussian priors are centred 
around the "true" parameter value, so a Gaussian prior imposing a variance of 
cr^ = 1 on a spin parameter x, with "true" value x = 0-95 will be a Gaussian 
prior around the value 0.95 and will allow a non-zero probability of finding a 
spin greater than 1. These are called normal true-parameter centered priors 
(NTC) in [23]. 

Slightly harder to treat, and not employed in the earlier works [7] and [8] 
are flat priors bounded by some region, that for example could correspond to 
positive mass conditions or maximal spin [23]. In this case care must be taken 
in choosing suitable parameters with which to impose flat priors. A flat prior 



on n, with upper cutoff iimax is different to a flat prior imposed on In(y^) with 
upper cutoff ln(/Xmax)- Also, flat priors in spin parameters /3 and a will not 
necessarily be the same as flat priors in xi and X2, when the relation between 
them is non-linear. 

The frequencies in (3) are the frequencies of the gravitational waves not 
the orbital frequencies of the massive objects. We will use the approximation 
that the gravitational waves have a frequency twice the orbital frequency of the 
objects. This usage is consistent with standard usage in the literature. In the 
following we will follow standard practice and mainly take /. — lOHz and /, as 
twice the frequency of a test-mass orbiting a Schwarzschild black hole of mass 
equal to the total mass of our system at the Innermost Stable Circular Orbit 
(ISCO), that is /^ = l/{6^/'^TrM). 

In order to solve (7) we take note of the fact that the integral, /, of a multi- 
dimensional Gaussian, with respect to a particular parameter 0p gives error 
functions 



,-r,,Ae,Ae,/2 f ^-r„Ael/2-rp,A0pA9,^^^Q^ 



= ^e^i^^^,-Tp^^,,/^„)A9.Ae,/2 (Erf[i^(A0;""=^)] - Erf[F(A6i™")]()l,0) 

where /, J are index labels containing the particular value p and i,j are index 
labels that do not. The argument of the error function in this case is given by 

When these integration limits, A6'™°^ and A^™*", are taken to ±cx) the differ- 
ence in the error functions just evaluates to two and the overall effect of the 
integration is to project out the direction Tpi in the Fisher metric F^ , up to an 
overall factor which cancels out in the posterior [24] . 

Multiple parameters can be marginalised in this way. Projecting out all 
parameters except one leaves just the variance that can be obtained directly 
by inverting the Fisher matrix. The effect of a flat prior will effectively be to 
change the limits of integration from ±oo to whatever maximum and minimum 
value the prior allows. In this way flat priors can be included, although multiple 
flat priors cannot be integrated analytically [23]. 

3 2.5PN Spinning Waveform 

The induced strain amplitude on the interferometer over time is related to the 
strain amplitude of the gravitational wave h(t). In order to compute Fisher ma- 
trix elements we need the Fourier transform of h{t). The Fourier transform is 



most efRciently performed analytically using the Stationary Phase Approxima- 
tion (SPA), which will be valid when the amplitude is changing slowly relative 
to the frequency of the wave [7]. The waveform in the Stationary Phase Ap- 
proximation is 

/i = A/-^/6e^*. (12) 

Schematically the phase can be written as 

vp, V]/„ Xj/o v]/. 

* = 2nfU + *o + -^ + -^ + -^ + -^ + *5 log(/), (13) 

where M is the total mass. The new term appearing at 2.5PN order is \l/5, 
although there are also new terms appearing in V^o and differences between the 
"^4 computed by Arun et al. [10] and the one used by Poisson and Will [8] 
due to self-spin interaction effects for the two spinning bodies. These self-spin 
corrections to the 2PN phase were first computed in [12]. The form of these 
expansion coefficients, as given in Arun et al. [10] are; 

.o^-..-M(f.Iog(^)), (14) 

*^ - T28(dlF^' ^''^ 

^,^ 3715 55m^/^ 

' " 327r2/3^3/2 8^^3/2 ' ^ ' 

15A^^/3 / 3058673 5429 / ^ yVz 617 / m V"^ \ 
'* ~ 647ri/3^2 1^1016064 ^ 1008 \m) ^ 144 IaI / "V' ^ ^ 

_^{My^^ f3m45n 657r^_^^^5/2 
where, 
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(X1+X2), (20) 
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V 4536 162 \m) ^ Ts" \m) 



(X1+X2) 
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V 4536 18 \M/ J \ \MJ 

_ M1M2 
^^ M1 + M2 



1-4(t7) (X1-X2), (22) 

(23) 



^ (M,M.)3A^ 

(Mi+M2)i/5- ^^^^ 

These are components of a Taylor F2 waveform in the classification of [25] . 
V^o contains both the phase at coalescence, the term — 7r/4 coming from the SPA 
and spin terms coming from 2.5PN order. Spin terms also occur in the last three 
terms ^3, ^4 and ^5. ^5 only contains spin-orbit contributions in the waveform 
of [10]. With seven terms with different functional dependencies on / one might 
expect to be able to reasonably determine seven parameters (including the time 
at coalescence coming from the linear / term.) To what extent this expectation 
is borne out depends on the details of the functional dependence of the chirp 
times and on the precise / dependence through the inner product (3). If the 
parameters are chosen such that the waveform depends linearly on them then 
the resulting Fisher matrix elements will be independent of the true parameter 
values. This however, is not the standard practice which we follow here. 

Cutler and Flanagan [7] considered only spin-orbit contributions at 1.5PN 
order. Poisson and Will [8] considered spin-orbit and spinl-spin2 contributions 
at 2PN order, only the first term of (21). In addition, using the waveform of 
Arun et al. [10] we consider spinl-spinl, spin2-spin2 terms at 2PN and spin- 
orbit terms at 2.5PN. In the following, unless stated otherwise, "2.5PN" refers 
to the waveform of [10] and "2PN" refers to the waveform used in [8]. 

To get some impression of the relative weight of the different expansion 
coefficients we present in Tables 1-5 their numerical values for a range of different 
parameters. We use units where the mass is given in seconds, so that Mq = 
4.93 X 10~^ seconds and frequencies are measured in Hertz. The coefficients 
have been rounded to three significant figures. The number of cycles refers to 
the total number of cycles the gravitational wave goes through between /. and 
/^. As expected the values of ^1, ^2 and ^3 are the same for both the 2.5PN 
and 2PN cases. In all cases the values of ^1, ^2 and ^3 are the same as the 
2PN versions and ^4 takes the same value for zero spin. There is a difference in 
the total number of cycles of the gravitational wave of up to 10 cycles. In some 
cases, for fairly large frequencies, / ^ lOOHz, the ^5 log(/) term gives a larger 
contribution to the phase than lower terms such as ^i/~^'^. 



4 Sensitivity curves 

The expected sensitivity of uncompleted interferometers can be simulated. These 
simulated sensitivites are typically obtained by modelling the thermal noise 





fo 


*i 


1'2 


*3 


*4 


*5 


no. cycles 


2PN 


-0.785 


66200 


1950 


-1030 


64.1 





227.6 


2.5PN 


-104 


66200 


1950 


-1030 


64.1 


14.5 


220.4 



Table 1: Comparison of 2PN and 2.5PN waveform coefficients for the parameter 
values Ml = M2 = IOMq, xi = X2 = and </-c = [f^sco = 220Hz). 
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Table 2: Comparison of 2PN and 2.5PN waveform coefficients for the parameter 
values Ml = 15Mo, M2 = IOMq and Xi = X2 = (fjsco = 176Hz). 
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Table 3: Comparison of 2PN and 2.5PN waveform coefficients for the parameter 
values Ml = I5M0, M2 = lOM© and xi = 0.95, X2 = 0.9 {f,sco = 176Hz). 
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Table 4: Comparison of 2PN and 2.5PN waveform coefficients for the parameter 
values Ml = 5Mq, M2 = 3M0 and xi = X2 = {f^sco = 550Hz). 
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1.39 
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54.6 


-0.27 


1157.5 



Table 5: Comparison of 2PN and 2.5PN waveform coefficients for the parameter 
values Ml = 5M0, M2 = 3M0 and xi = 0.95, X2 = if,sco = 550Hz). 
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which dominates at lower frequencies and photon shot noise which dominates 
at higher frequencies. Seismic noise can be included by adding a sharp cutoff 
at a low frequency. The initial LIGO interferometers were able to attain their 
design sensitivity over large parts of the frequency range [27]. With the same 
waveform and a fixed signal to noise ratio, the shape of the sensitivity curve 
determines the parameter errors. The power spectral density (PSD) used by 
Cutler and Flanagan [7], with /o = 70Hz, is 



Shif) = 3 X 10- 



-48 



2 1 




(25) 



The sensitivity curve used here for advanced LIGO [28], with /o = 215Hz, is 



Shif) = 10- 



JO/ V /o , 



111 



0.5 



iiT 



l + 0.5(^ 



(26) 

The curve used by Vallisneri [23], which is the initial LIGO curve from Table 
IV of [25], with /o = 150Hz, is 



Shif) = 9 X 10 



-46 



4.4., I 1))-% 0.16 (i.)"" + 0.52-^0.32 a 



(27) 



A projected curve for the advanced Virgo detector can be obtained from [26] 



5 Results 

We choose to adopt parameters for the waveform tc, 4>c, M, A^, Xi ^^id X2, with 
tc measured in milliseconds. The choice of these parameters is largely dictated 
by simplicity and convention. Other combinations of these parameters, such as 
the symmetric mass ratio, 77, or the spin functions /3 and a, as used in [8], can 
be solved for using the propagation of errors formula [21]. The linear spin terms 
/3, a and 7 are not independent parameters at 2.5PN order and for our purposes 
it is easier to use xi and X2 directly. 

The waveforms do not distinguish the two black holes if their masses and 
spins are equal. If we take parameters to be (tc, '?!>c,ln/u,ln7\/J,xi,X2) then the 
Fisher matrix inversion encounters problems when Mi = M2 and Xi = X2 
as the two objects cannot be distinguished and the parameters xi ^^d X2 are 
degenerate leading to a Fisher matrix with zero determinant, so we deliberately 
avoid this limit by taking slightly different mass values. 

Using the same parameters as [8] with /? and a and their 2PN waveform 
we are able to reproduce their results up to ~ 3% accuracy. It is not clear 
what this discrepancy is due to. The numerical integration has been checked 
for convergence and reproduces successfully the results of [7] using the 1.5PN 
waveform. The main source of error in our results in the numerical inversion of 
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At, 


A0C 


Am/m 


AM/M 


Axi 


Ax2 


2.5PN 


37.0 


191.5 


7.74 


0.0975 


261.1 


430.4 


2PN 


22.4 


248.3 


21.1 


0.157 


14.0 


44.5 



Table 6: Comparison of parameter errors for 2PN and 2.5PN waveforms with 
parameter values Mi = 15Mq, Ma = IOMq, xi = 0.95 and X2 = 0.9. 



the Fisher matrix, although for most of our results this is accurate to one part 
in 10^. We are also able to reproduce the 2PN results of [23] to the accuracy 
reported there using the sensitivity curve (27) and a lower frequency cutoff at 
f. = 40Hz and the results for the 2.5PN non-spinning waveforms of [11] for a 
lower cutoff at /. = 20Hz. 



5.1 Differences between 2PN and 2.5PN 

We begin by comparing the 2.5PN waveform of [10] with its seven expansion 
coefficients and the 2PN waveform used in [8] with only six expansion coefficients 
and slightly different functional dependence in ^o and ^4. In the following, we 
will denote the waveform used in [8] by "2PN" and the symbol A9 is used 
to denote the root mean squared error relative to the mean - the standard 
deviation. This is not the same as the difference to the mean used in Eqn. (4), 
but is consistent with standard notation in the literature. 

Parameter error values for true parameters Mi = 15M0, M2 = IOMq, 
Xi = 0.95 and X2 = 0-9, without any priors are given in Table 6. The error on 
H is considerably better for the 2.5PN waveform with a good improvement in 
M too. The worsening in the spin parameters is extreme, although neither the 
2.5PN nor 2PN waveforms are close to the physically expected bound \x\ < 1. 

For the same true parameter values with Gaussian priors on both the spins 
a^i = 1 and ax2 = I7 the error values are given in Table 7. There is still an 
approximately factor three improvement in In^, but a loss of accuracy in In(A^). 

In Table 8 we give error values for example parameter values of a light, 
rapidly spinning black hole and non-spinning neutron star system, with Mi = 
5M0, M2 = I.4M0, xi = 0.95 and X2 = without priors. The errors are once 
again better at 2.5PN for the mass parameters and worse for the xi spin pa- 
rameter, although this time better for X2- The error difference is still significant 
in fi even for these low masses. 

Errors for the same true parameters with Gaussian priors on both the spins 
a^-^ = 1 and a^^ — 1 are given in Table 9. There is still some difference in the 
parameter estimation errors. From this we conclude that using 2.5PN waveforms 



12 





At, 


A(/), 


Afi/fi 


AM/M 


Axi 


Ax2 


2.5PN 


7.44 


13.6 


0.61 


0.0050 


0.77 


0.92 


2PN 


1.66 


16.8 


2.06 


0.0019 


0.98 


1.00 



Table 7: Comparison of parameter errors for 2PN and 2.5PN waveforms with 
parameter values Mi — 15Mq, M2 ~ IOMq, xi = 0.95 and X2 = 0.9 and 
Gaussian priors on both the spins a^i = 1 and a^^ = 1. 





At, 


Acj), 


Afi/fi 


AM/M 


Axi 


Ax2 


2.5PN 


3.70 


23.3 


0.69 


0.0039 


1.4 


5.0 


2PN 


2.84 


54.2 


1.04 


0.0048 


0.85 


10.4 



Table 8: Comparison of parameter errors for 2PN and 2.5PN waveforms with 
parameter values Mi = 5Mq, M2 = 1.4Mq, xi = 0.95 and X2 = 





At, 


Acj), 


Afi/fi 


AM/M 


Axi 


Ax2 


2.5PN 


3.07 


19.1 


0.54 


0.0028 


0.60 


0.96 


2PN 


0.79 


5.42 


0.22 


0.0016 


0.18 


0.99 



Table 9: Comparison of parameter errors for 2PN and 2.5PN waveforms with 
parameter values Mi = 5M0, M2 = I.4M0, xi — 0.95 and X2 = and Gaussian 
priors on both the spins a^^ = 1 and a^^ — 1- 
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At, 


A0C 


Am/m 


AM/M 


Axi 


Ax2 


2.5PN 


37.0 


191.5 


7.74 


0.098 


261.1 


430.4 


2PN with self-spin 


22.4 


763.8 


21.1 


0.157 


502.4 


841.0 


2PN 


22.4 


248.3 


21.1 


0.157 


14.0 


44.5 



Table 10: Comparison of parameter values for parameters Mi — 15Mq, M2 = 
IOM0, xi = 0.95 and X2 = 0.9 and ^5 = 0. For clarity, "2PN" should be here 
understood as the 2PN waveform used in [8]. "2PN with self-spin" is 2PN with 
the spin-orbit and spin-spin terms considered in [8] and also the self-spin terms 
of [10], which is equivalent to the 2.5PN waveform with \I'5 = 0. 



of [10] will give noticeably different results to using just the 2PN waveform of 
[8]. It is worth remembering that these results are based on the assumption 
that the true waveform is a 2.5PN waveform for the 2.5PN numbers and a 2PN 
waveform for the 2PN ones, neither of which is likely to be true of the full 
General Relativistic waveform. 

To see what causes the large errors differences obtained from the 2PN and 
2.5PN waveforms used here, we can artificially remove the log(7rM/) dependent 
term in the 2.5PN waveform while keeping the self-spin modifications to ^4. 
The resulting errors for a rapidly spinning black hole binary system with true 
parameters Mi = I5M0, M2 = IOMq, xi = 0.95 and X2 = 0.9 are displayed 
in Table 10. The errors in the xi &n.d X2 parameters are much larger, larger 
even than at the full 2.5PN level and there is also a corresponding increase in 
the error in (pc- This suggests that much of the difference between the 2.5PN 
waveform of [10] and the 2PN of [8] is due to the self-spin interaction terms at 
2PN not included in [8]. 

The larger error ranges noted for the spins may also be partly due to the xi 
and X2 dependent terms that appear in ^I^o at 2.5PN level and cause confusion 
with the (j)c term. If wc use ^0 as a parameter, rather than (j), we obtain results 
displayed in Table 11. The errors for parameters other than (Jjc/Psig are seen 
to be unchanged. We conclude that mass and spin parameter errors cannot be 
changed by using ^0 as a parameter instead of (pc- This is consistent with the 
findings in [11]. 

There is still a large difference from the results obtained with the 2PN wave- 
form of [8] . The Fisher matrix is more degenerate in the spins when the self-spin 
couplings are included, meaning it is harder to distinguish the spin values. This 
large difference in the errors may be partly due to the technical fact that the 
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At, 


A*o 


Am/m 


AM/M 


Axi 


Ax2 


full 2.5PN 


37.0 


1202.4 


7.74 


0.098 


261.1 


430.5 


2PN with self-spin 


22.4 


248.3 


21.1 


0.157 


502.4 


841.0 



Table 11: Comparison of errors for parameters Mi — 15Mq, M2 = IOMq, 
Xi = 0.95 and X2 = 0.9 and using vJ/q as a parameter, instead of (j),- 



extra self-spin terms are quadratic in the spin parameters and the Fisher for- 
malism works best with the Linear Signal Approximation [23]. However, the 
second derivatives of the waveform phase with respect to the spins was found 
to be small with respect to the first order terms (by a factor ^ 0.2) and the 
expansion terminates at second order. The waveform phase is also not linear in 
the chirp mass or symmetric mass ratio, even at the 2PN order considered in 
[8]. 

The effect of using large, near extremal true spin values can be compared to 
assuming their true value is zero. For Mi — 15Mq, M2 — IOMq, Xi = X2 = 
we find the the errors in In/n and InTM increase to 17.13 and 0.14 respectively, 
but the errors in xi smd X2 are only 59.1 and 49.7. This is much lower than the 
261.1 and 430.4 found when they are assumed to be nearly maximally spinning. 
Even when the true values of the spins are zero, including the errors due to the 
spins has a large detrimental effect on our ability to measure the masses exactly. 
In the case of Mi — 15Mq, M2 = IOMq, Xi = X2 = 0, if we assume that the 
spins are exactly zero, and do not marginalise over them then the measured 
errors in In/x and In A^ are only 0.020 and 0.0031, an improvement of roughly 
three, respectively two, orders of magnitude. However, this increase in accuracy 
will require outside information from other observations or prior restrictions and 
cannot be justified from the gravitational waveform alone. 

5.2 Effect of frequency ranges and priors 

The final state of the merger of two compact objects is likely to be a spinning 
black hole. For progenitors with large, aligned spins, the final black hole is 
likely to have a large spin itself. In this case it is questionable to what extent 
the ISCO of a non-spinning Schwarzschild black hole is a suitable indication of 
where the FN inspiral approximation becomes unreliable. A near-extremal Kerr 
black hole has an ISCO at r ^ M, as opposed to r ^ 6M for a non-spinning 
hole. A spinning black hole can support test particles in quasi-circular obits at 
much higher frequencies than a non-spinning black hole of the same mass. The 
asymptotic time period of a test particle in a circular orbit in the equatorial 



15 



plane of a Kerr black hole is T = (27r(a + \Jr'^ /M)) and hence the gravitational 
wave frequency at extremality is approximately fj — l/(27rM). 

There is still a lot of uncertainty about where PN approximations break 
down and the Schwarzschild ISCO has traditionally been used as a conservative 
approximation. It is possible [18] that the PN waveform for a rapidly spinning 
system is valid to a significantly higher frequency than this. We can look at the 
approximate effect of increasing the final frequency by allowing fi to increase by 
a factor of ten. The frequency difference between the true Schwarzschild ISCO 
and the true extremal Kerr ISCO is closer to a factor of 7, but as there is still 
uncertainty in exactly where the true limit should be applied, we adopt here a 
factor of ten for simplicity. 

Prior information can also reduce the uncertainty in parameters by a sig- 
nificant amount. In [7] and [8] Gaussian priors were imposed to simulate in 
a simple analytical fashion the effect of exact physical priors on the spins of 
compact objects. To identify the effect of different priors and changing the in- 
tegration range, in Table 12, with true parameters Mi = SMq, M2 = SM©, 
Xi = 0.95, X2 = 0.9 and the SNR normalised to 10, we implement several differ- 
ent choices of Gaussian spin priors cr^^ ^ — 1 ^^d a Gaussian coalescence phase 
prior CT0^ = TT. A longer inspiral waveform with SNR normaised to 10, will have 
a lower amplitude than a shorter waveform in the corresponding overlap region. 

An increase in f^ by a factor 10 such as might occur in the proposal of [18], 
can lead to a substantial gain in accuracy for unpriored parameters, although 
the effect diminishes when priors are applied. The Gaussian priors a^^ = tt 
and (T^j 2 ~ 1 contain nearly all the information about these parameters for 
integrating both to the Schwarzschild ISCO and ten times the Schwarzschild 
ISCO. The spin priors have a greater effect on both Aln/i and AlnTM than the 
(j)c but the combination of both improves the mass estimation by about a factor 
twenty. A prior on 0c was rejected in [23] since the coalescence phase, (j)c, is 
able to absorb the phase shifts due to changes in the other parameters. 

Another prior one could impose is a prior on the masses. The maximum 
possible value of the symmetric mass ratio rj = /i/(Mi -I- M2) is 0.25 which 
occurs for equal masses. This constraint is simply related to the requirement 
that the inspiralling objects have real positive masses. It can be translated into 
a maximum value on fj. of M/4 for a fixed total mass M . Such a prior can be 
implemented as a flat prior in the allowed range using the above mentioned like- 
lihood function. The lower bound on the reduced mass at /i = is automatically 
implemented by using the parameter In/z with no bound below. 

Imposing the upper bound on In/i it was found that the variance on the 
Xi parameter fell from 261.1 to only 260.5 for the parameters Mi ~ I5M0, 
M2 = IOM0, xi — 0.95 and X2 = 0.9. This is still far from giving a meaningful 
constraint and cannot exclude the unphysical regime xi > 1- However, the 
variance in this case does not give equal probabilities on either side of the 
mode, but is heavily skewed towards lower values and the mode itself is shifted 
away from the "true" value. A flat prior in In/i corresponds to a logarithmic 
prior in the parameter fi since it implies that each decade of fi has equal prior 
probability. 
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At, 


A0e 


A/i/M 


AM/M 


Axi 


Ax2 


no priors 


4.9 


13.0 


1.31 


0.0076 


17.2 


31.1 


no priors 


1.6 


9.3 


0.64 


0.0043 


10.9 


20.0 


spin priors 


2.6 


12.4 


0.35 


0.00094 


0.55 


0.89 


spin priors 


1.1 


5.7 


0.18 


0.00050 


0.50 


0.88 


(f>c prior 


4.2 


3.1 


1.26 


0.0076 


17.2 


31.1 


spin and (f)c prior 


0.68 


3.0 


0.10 


0.00044 


0.49 


0.88 


spin and (fie prior 


0.56 


2.8 


0.099 


0.00040 


0.49 


0.87 



Table 12: Effect of various combinations of the Gaussian priors a^-^ 2 — 1 ^^'^ 
a^^ = TT, and increasing the final frequencies for the parameters Mi = 5Mq, 
M2 = 3Mq, xi = 0.95, X2 = 0.9 and SNR normahsed to 10. 
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A more impressive constraint can be obtained by using the parameter 77 
directly instead of In/x. In this case the cutoff at r; = 0.25 can be imposed 
directly as a flat prior in rj between and 0.25. For high masses Mi = 15Mq, 
M2 = IOMq, and near extremal spins xi — 0-95 and X2 = 0.9, Axi is reduced 
from 261.1 to 27.6 and Ax2 is reduced from 430.4 to 47.3 (There is also a small 
shift in the Fisher matrix values, presumably due to the non-linear change of 
parameters from In/i to rj.) For lower masses. Mi = 5Mq, M2 = SMq, with 
true spins set to zero, the reduction in Axi is from 22.9 to 6.9 and the reduction 
in Ax2 is from 34.9 to 12.5. For the unequal mass BH-NS case. Mi = IOMq, 
M2 = 1.4Mo with xi = 0.95 and X2 = 0, Axi reduces from 1.9 to 1.2 and Ax2 
from 15.5 to 10.6. These reductions are still not sufficient to place the spins 
inside the physical bound and there does not seem to be any strong reason for 
preferring a flat prior in 77 to a flat prior in In fi other than the more direct 
relation to the cutoff at rj = 0.25. 

Although it does not seem possible with the 2.5PN waveform to assign well 
measured spin values to the individual components of the binary, one of the 
main reasons for this is that the spin parameters xi and X2 are not very well 
distinguished by the 2.5PN aligned spins waveform. The two spin parameters 
are highly degenerate. However, the error ellipse given by the likelihood function 
obtained by marginalising over all parameters except the spins is often very thin 
in one direction. This means that there is a particular linear combination of 
the spin parameters that is much better constrained by the waveform than the 
individual spins. In [29] it was suggested that this combination might be 

X-^-Xi + ^-X2, (28) 

where S = (Mi - M2)/(Mi + M2). In the case with Mi = 15Mo, M2 = IOMq 
we have S — 0.2 and thus x = 0.6xi + 0.4x2- 

With a two-dimensional likelihood function, obtained by marginalising over 
the other parameters, the preferred direction can be found by finding the eigen- 
vectors of the associated Hessian matrix. Numerically with Mi = 15Mq, 
M2 = IOM0, xi = 0.95 and X2 = 0.9 we find x = 0.86xi + 0.52x2- The 
error on this linear combination is Ax = 0.91 and the maximum value it can 
attain, when both spins are at their maximal physical value, xi = X2 ~ Ij is 
1.37. For the same high masses, but xi = X2 = we find x = 0.60xi + 0.80x2. 
The error on this combination is Ax — 33.2 and the maximum value it can 
attain is 1.4. 

For lower masses, such as Mi — 5, M2 — 3 with near maximal, aligned 
spins, xi — 0.95 and X2 — 0.9, we find x — 0.88xi + 0.48x2 as the optimally 
measured linear combination of spin parameters. The error on this combination 
is just Ax = 0.31 and its maximum possible physical value is 1.36. With no 
true spin we obtain instead x = 0.84xi + 0.55x2 as the optimally measured 
linear combination. This combination has a measurement error of Ax = 1.28. 
The maximum value that x can attain if both spins are at their maximal phys- 
ical value is 1.38. This suggests that in certain parameter ranges, particularly 
lower masses, such as neutron stars or light black holes, it may be possible to 
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At, 


A'^c 


Am/m 


AM/X 


Axi 


Ax2 


6 params 
X2 unknown 


4.7 


28.7 


1.52 


0.0082 


21.4 


30.6 


5 params 
X2 known 


2.9 


19.6 


0.50 


0.0016 


1.85 


0.0 (exactly) 



Table 13: Comparison of estimation of six parameters versus five parameters 
for true parameter values M\ = 5, M2 = 3, xi = 0.95 and xi = 0. 





At, 


A0e 


Am/m 


AM/M 


Axi 


Ax2 


6 params 
X2 unknown 


7.0 


15.0 


0.84 


0.0079 


1.9 


15.5 


5 params 
X2 known 


4.7 


13.2 


0.38 


0.0025 


0.28 


0.0 (exactly) 



Table 14: Comparison of estimation of six parameters versus five parameters 
for true parameter values Mi — 10, M2 — 1.4, xi — 0.95 and X2 = 0. 



demonstrate that the system contains a non-zero spin within a physical bound, 
although separate spins cannot easily be assigned to the individual objects. 

In the case that extra external information is able to precisely fix the spins 
of one of the objects, perhaps the lighter object, then the parameter estimation 
problem becomes essentially five dimensional. This might happen if one of the 
objects is a pulsar, or a neutron star with known equation of state that prevents 
it from having any meaningful spin angular momentum. In this case the spin 
of the heavier object can be measured much more accurately without the need 
for priors (although the spin of the other object is essentially a delta function 
prior). The low mass case Mi = 5, M2 = 3 with true parameters xi — 0.95 
and X2 = are compared in Table 13 and a black hole-neutron star example, 
with Ml = 10, M2 = 1.4 and true parameter xi — 0.95 and X2 = in Table 
14. In the case of a low mass companion, with the lighter spin known exactly, 
it may be possible to determine that the spin of the heavier object lies within 
the physical range to Ict, without the use of any priors. 

We also present a comparison of the 2.5PN inspiral waveform and the phe- 
nom waveform (which includes 2.5PN spin terms) from [30] in Table 15. This is 
necessarily only a four-dimensional parameter set, as the way spins arc treated 
is slightly different in the two cases. Effectively here we are assuming that the 
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Aic 


A(/), 


Am/m 


AX/X 


Ml = 15Mo, M2 = IOMq 

2.5PN waveform 

Xi = X2 = 


1.8 


3.04 


0.020 


0.0031 


Ml = 15M0, M2 = IOM0 
phenom waveform 
X = 0, /, = ISCO 


2.3 


3.69 


0.027 


0.0039 


Ml = 15M0, M2 = IOMq 

phenom waveform 

X = 0, /, = 0.15/M 


1.6 


2.83 


0.022 


0.0035 


Ml = 5M0, M2 = 3Mq 

2.5PN waveform 

Xi = X2 = 


0.56 


1.60 


0.0070 


0.00037 


Ml = 5Mq, M2 = 3M0 
phenom waveform 
X = 0, /, = ISCO 


0.51 


1.34 


0.0059 


0.00036 


Ml = 5M0, M2 = 3Mq 

phenom waveform 

X = 0,/,= 0.15/M 


0.41 


1.15 


0.0052 


0.00033 



Table 15: Comparison of PN inspiral waveform and phenom waveform of [30] 
for various true parameters and integration ranges. 



true spins are zero and known to be exactly zero. 

With the merger and ringdown modeled we can extend the integration 
range weh beyond the ISCO frequency, here we can take it to 0.15/M (sug- 
gested by Frank Ohme) which corresponds to 1218Hz instead of 176Hz for the 
Schwarzschild ISCO with total mass 25M0 and 3800Hz instead of 550Hz for 
the total mass 8M0 case. When extending the frequency cutoff, to ensure a fair 
comparison, we can fix the amplitude and allow the SNR to increase. In the 
high mass example, there is some increase in the SNR from 10 to 11.72 but a 
slight deterioration in the mass parameters' accuracy. In the low mass example 
the SNR only increases from 10 to 10.0147 but there is an increase in error on 
In/i and no meaningful change in the chirp mass. 

This analysis has ignored the spins, effectively assuming that they are known 
and zero, but it does not appear that there is much to gain from using this 
particular phenom waveform, even after the frequency range is greatly extended 
and in this way the SNR somewhat increased for the same amplitude signal. 
This interesting area though deserves more study as better phenom and hybrid 
Post-Newtonian and numerical waveforms become available. 
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At, 


Acj), 


Afi/fi 


AM/M 


Axi 


Ax2 


SfadvLIGO 
no priors 


4.7 


28.7 


1.5 


0.0082 


21.4 


30.6 


SfadvLIGO 
spin priors 


1.63 


10.6 


0.26 


0.00086 


0.89 


0.95 


SffnitLIGO 
no priors 


12.2 


38.7 


8.5 


0.074 


184.6 


294.1 


SffnitLIGO 

spin priors 


1.5 


10.5 


0.28 


0.0019 


0.97 


0.99 


SfCutFlan 
no prior 


6.9 


42.7 


1.45 


0.0066 


14.5 


18.6 


SfCutFlan 
spin priors 


2.3 


12.8 


0.28 


0.0010 


0.95 


0.97 


SfadvVirgo 
no prior 


3.6 


23.1 


1.16 


0.0061 


15.7 


22.2 


SfadvVirgo 
spin priors 


1.4 


9.9 


0.24 


0.00079 


0.86 


0.93 



Table 16: Comparison of different sensitivity curves both with and without 
Gaussian spin priors for true parameters Mi = 5Mq, M2 = 3Mq, xi = 0.95, 
X2 = 0. 



5.3 Shape of sensitivity curve 

It is also interesting to compare the effect of different shapes of sensitivity curve. 
The sensitivity curves used in [8], [23] and here are all different, see Eqns (25- 
27). Using the low mass parameters Mi — 5Mq, M2 = 3Mq, xi — 0.95, 
X2 = 0, with a normalised SNR= 10, an initial frequency f- = lOHz and /^ at 
the Schwarzschild ISCO we obtain the results displayed in Table 16. 

We see that the advanced Virgo curve appears to be the best both with and 
without priors. This is likely due to the wider, but flatter shape of the advanced 
Virgo sensitivity curve. These sensitivities refer to singular detectors operating 
in isolation. When operating together with combined analyses the results will 
be better, partly through an improvement in effective SNR [7]. 

Finally, we list some results from various different regions of the parameter 
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space in Table 17. With both spins anti-ahgned to the orbital angular momen- 
tum the spin uncertainty becomes very large. 

6 Conclusions 

We have obtained a number of results related to parameter estimation using 
Post-Newtonian 2.5 order spinning waveforms. The most important result is 
that measurement errors differ significantly from previous results obtain by Pois- 
son and Will [8] at 2PN order. A large part of this difference is attributable to 
the self-spin interaction terms at 2PN order not included in the analysis of [8] 
or in analyses for space-based missions [20]. However, it also strongly suggests 
that 2PN spinning waveforms are not sufficient for reliable parameter estimation 
from spinning systems and possibly 2.5 PN order is not sufficient either. Using 
the best available 2.5PN waveform leads to a Fisher matrix that is more degen- 
erate in the spins than previous studies. This also means that small changes in 
the waveform can have large effects on the computed errors. Source multipole 
moments needed to obtain the spin contributions for gravitational waves have 
recently been computed at 3PN order [31] and it is hoped that the results re- 
ported here will give extra impetus to putting these contributions into a "ready 
to use" form necessary for data analysis and parameter estimation. Higher order 
spinning waveforms may also be needed beyond 3PN. 

We have also explored a number of issues with the 2.5PN waveform. Implic- 
itly the analysis presented here assumes that nature corresponds exactly to this 
2.5 PN order. The model of an exactly aligned, 2.5PN order, spinning wave- 
form of this form is a prior assumption in the analysis. Furthermore the analysis 
assumes that an individual detector will have a stationary Gaussian noise spec- 
trum given, for example, by the modeled approximation (26) and ignores any 
other systematic errors in the observational and data manipulation processes. 
In light of these facts, and the belief that nature is unlikely to be exactly 2.5 
PN, it is difficult to give precise estimates for the parameter errors that can 
be obtained from real advanced second generation interferometers. The values 
given above should taken as purely indicative of how signals might behave in 
certain parts of parameter space and which parts of parameter space might be 
amenable to constraining which parameters. Fisher matrix calculations have 
been compared with Monte Carlo methods in [19] and there it was found that 
Fisher matrix techniques can in cases underestimate the errors by a factor of 
two. 

We have shown that while the mass parameters are often well constrained 
using the 2.5PN waveform, allowing in cases black hole candidates to be dis- 
tinguished from neutron star candidates, the spin information is only weakly 
constrained. Only in the case where one of the spins is known exactly from 
external information can the spin of the other object be constrained within the 
physical bound |xi| < 1. 

The detectors have limited sensitivity to intrinsic spin at 2.5PN level and 
information about the spins is still likely to be dominated by prior information. 
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At, 


A(j)c 


Am/m 


AM/M 


Axi 


Ax2 


5Mo, 3Mq 

XI = 0, X2 - 

no priors 


4.3 


80.0 


1.80 


0.0089 


22.9 


34.9 


5M0, 3Mq 

Xi =0.95, X2 =0.9 

no priors 


4.9 


13.0 


1.3 


0.0076 


17.2 


31.1 


5Mo, 3Mq 

XI = 0.95, X2 = -0.9 

no priors 


4.3 


55.1 


1.73 


0.0087 


24.9 


27.0 


5Mo, 3Mq 

XI - 0, X2 = 

spin priors 


0.6 


9.6 


0.35 


0.0024 


0.86 


0.94 


15M0, IOMq 

Xi=0.95, X2=0.9 

no priors 


37.0 


191.5 


7.7 


0.097 


261.1 


430.4 


15M0, IOMq 

Xi = 0, X2 = 

no priors 


27.0 


536.6 


17.1 


0.14 


59.1 


49.7 


15M0, IOMq 

XI - 0.95, X2 = -0.9 

no priors 


26.8 


373.8 


17.3 


0.14 


607.9 


745.6 


15M0, IOM0 

XI = -0.95, X2 = 0.9 

no priors 


18.9 


1160 


24.9 


0.17 


991.4 


1051.7 


15M0, IOM0 

Xi = -0.95, X2 = -0.9 

no priors 


21.1 


1623.7 


23.1 


0.17 


994.0 


1398.4 


15M0, IOM0 

XI = 0.95, X2 = -0.9 

spin priors 


2.2 


6.5 


0.13 


0.0033 


0.996 


0.997 


IOM0, IAMq 

Xi = 0.95, X2 = 

no priors 


7.0 


15.0 


0.84 


0.0079 


1.9 


15.5 


IOM0, IAMq 

XI = -0.95, X2 = 
no priors 


7.0 


149.0 


0.88 


0.0081 


11.4 


82.1 


IOM0, IAMq 

XI = 0, X2 = 
no priors 


6.5 


99.3 


1.0 


0.0088 


4.9 


35.0 



Table 17: Parameter errors for various different true parameter values for the 
2.5PN aligned spinning waveform and advanced LIGO sensitivity curve. Some 
values are also given with Gaussian priors on the spins of tr^j — 1 and tr^j = 1. 
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Certain information that is necessary in the parameter estimation is difficult to 
fold into the Fisher matrix formalism. More accurate prior information can be 
encoded into the likelihood using a variety of techniques [23] . The method used 
here employs a semi-analytical technique that greatly reduces computational 
complexity, but is limited to one flat prior and the exact choice of parameter 
is an important consideration in applying "flat" priors. The prior information 
that may be so important to obtaining reliable mass estimates, is still not well 
implemented. 

Assumptions about where the inspiral PN breaks down, using Kerr ISCO 
instead of Schwarzschild ISCO, does have some effect, but prior information is 
still dominant in these cases. Phenom waveforms do not have sufficient infor- 
mation or accuracy yet to completely alleviate this problem. The 2.5PN level 
expansion is unlikely to be sufficient for reliable extraction of spin parameters 
in General Relativity. As advanced hybrid waveforms are developed containing 
a full waveform from inspiral to ringdown in full General Relativity more work 
will be needed to study this interesting area. 

Several proposals for measuring deviations from General Relativity rely on 
measuring simultaneously the parameters needed to describe the binary using 
the waveform and also checking for deviations from the predictions of General 
Relativity. With non-spinning systems the waveforms are largely just a function 
of the two masses, which can be measured fairly accurately using only a few chirp 
times, leaving the remaining chirp times as a consistency check of the model 
predictions. For aligned spinning systems, there are four intrinsic parameters 
to be measured to describe the system and as we have seen the spins are not 
well determined, even within the physically expected range. It is possible that 
waveforms at 3PN and above will alleviate some of these problems encountered 
here and leave this to further work. 
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8 Appendix 

We give here for completeness examples of the covariance matrices, S*-'. (We 
thank Drew Keppel for the suggestion that this would be useful to include.) 
The square root of the diagonal elements corresponds to the measured accu- 
racy of each parameter and the correlation between parameters a and b can be 
calculated by the formula 

y^ab 

c , = , (29) 
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For the case of six parameters, with tc = 4>c = 0, Mi = 15Mq, M2 = IOMq, 
Xi = 0.95 and X2 = 0.9 wc have 

15.356 \ 
-82071.1 
-3310.16 

41.886 
-112369. 
185276. J 
(30) 
0, Ml = 15M0, M2 = IOM0 and the 
spins assumed to be exactly known at their values xi — 0.95 and X2 — 0.9, we 
have 



/ 0.00137 


-6.627 


-0.282 


0.00343 


-9.323 


-6.627 


36669.7 


1451.79 


-18.632 


49757.5 


-0.282 


1451.79 


59.831 


-0.744 


2008.48 


0.00343 


-18.632 


-0.744123 


0.00951 


-25.398 


-9.323 


49757.5 


2008.48 


-25.398 


68152.8 


l^ 15.356 


-82071.1 


-3310.16 


41.886 


-112369 



^6 



With only four parameters, tc 



^4 



/ 5.906 X 10-6 
0.0102 
-0.000107 
\ 6.272 X 10-6 



0.0102 
18.091 
-0.191 
0.0115 



-0.000107 
-0.191 
0.00202 

-0.000124 



6.272 X 10-6 \ 
0.0115 
-0.00012 
9.141 X 10-6 J 



(31) 



The corresponding correlations are 



/ 1- 
-0.935 

-0.986 

0.950 

-0.965 

\ 0.964 



-0.935 
1. 

0.980 
-0.998 

0.995 
-0.996 ^ 

1. 

0.992 

-0.981 

0.854 



V 



-0.986 

0.980 

1. 
-0.987 

0.995 
-0.994 

0.992 

1. 

-0.997 

0.895 



0.950 
-0.998 
-0.987 

1. 
-0.998 

0.998 



-0.965 
0.995 
0.995 

-0.998 
1. 

-0.999 



0.964 \ 
-0.996 
-0.994 

0.998 
-0.999 

1- / 



-0.981 0.854 
-0.997 0.895 

1. -0.915 
-0.915 1. 



The parameters are all highly correlated, especially the spins. 



(32) 



(33) 
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